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Abstract 



We propose an evolutionary mechanism of phyllotaxis, regular arrangement of leaves on a plant stem. It is shown 
that the phyllotactic pattern with the Fibonacci sequence has a selective advantage, for it involves the least number of 
£N| phyllotactic transitions during plant growth. 
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1. Introduction 

Phyllotaxis, regular arrangement of leaves on a plant 
stem, has since long attracted the minds of botanists, 
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mathematicians and physicists 
monly, alternate leaves along a twig execute a spiral with 
an angle of 1/2, 1/3, 2/5, 3/8 of a full rotation, or other- 
wise with a limit angle of 137.5 degrees (Fig. [TJ. Patterns 
with the other fractions are also observed, though uncom- 
monly. Hence phyllotaxis is regarded not as a universal 
'law but as a fascinatingly prevalent tendency ([J]). For 
decades, mathematical works have elucidated number the- 
oretical structure of phyllotaxis and deepened our under- 
standing of the subject significantly (0, m, M, M, H, 0, 
20, 14f). Nevertheless, there still remains a fundamental 
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problem of why and how only some specific numbers are 
favored by nature. This is a problem of natural science. 

Physical or chemical models studied thus far are mostly 
based on dynamical mechanism, according to which a phyl- 
lotactic pattern is supposed to be realized autonomously 
as an end result of its own mechanical or chemical dynam- 
ics 

(0, US HI Q !])• This is mechanical determinism. 

Mathematicians and physicists generally appreciate this 
point of view. On the contrary, biologists generally take 
the opposite viewpoint that there are so strong reasons 
for the plants to have genetic information that a diver- 
gence angle between adjacent leaves must be determined 
genetically. This is genetic determinism. In this case, we 
still have to make clear how the plants inherit the math- 
ematical characteristics, i.e., a static or statistic mecha- 
nism of phyllotaxis is asked for. When it comes down to 
it, there seems no hope to do without mathematics, for we 
have to explain why the divergence angle is programmed to 
take neither 130° nor 140° but the magic number 137.5°. 
Unfortunately, candidate explanations have been only de- 
scriptive and qualitative, or otherwise quantitative but 
teleol ogica l such that packing efficiency ([l7jp or unifor- 
mity ( jlBl HH, Hi]) is meant to be maximized resultingly. 
In any case, at this fundamental level, we have to regard 



phyllotaxis as a conundrum of theoretical physics (or the- 
oretical biology), for we have to deal with mathematics 
observed in the real world. The aim of this paper is to 
present a satisfactory static mechanism to be contrasted 
not only with the existing dynamical models but with the 
teleological explanations. 

We propose a new mechanism based on a growth model 
of a plant. The model is basically characterized by two pa- 
rameters, an initial (preset) value of the divergence angle, 
2nao, and a vertical range of repulsive interaction, n c . The 
model is based on the observation that a plant has inherent 
abilities (1) to arrange leaves primarily on a regular spiral 
with a constant angle of rotation 2jrao, and (2) to exert 
secondary torsions between the leaves within reach of n c . 
The former is consistent with experiments on apical meris- 
tems (d, H^l). The latter is supposed to operate in a vas- 
cular system ((HI). The secondary interaction reinforces 
the regularity of the helical arrangement, thereby the ob- 
served angle a generally undergoes changes from the initial 
value ao. With the vertical range n c regarded as a growth 
index, it is shown that a goes through stepwise transitions 
between phyllotactic fractions (PF) in the course of the 
growth (Fig. In fact, many plants progress through 
distinct phyllotactic transition (PT) during early stages of 
development. It should be rather remarked that the sec- 
ondary changes, essential for us, are usually noticed but 
often disregarded as irrelevant (as secondary). To outline 
the proposed mechanism, consider a population of ran- 
dom samples with all possible values of the preset angle 
ao, and let them grow according to the model. As they 
grow, each sample will exhibit its own PTs respectively, 
depending on a.Q. With all the grown samples, we leave 
from formal analysis. At this point, we appeal to biologi- 
cal reasoning according to natural selection, by which all 
but that with the golden angle 2irao = 137.5° turns out to 
be eliminated. The selected sample is favored in nature, 
for it is structurally the most stable because it undergoes 
the least PTs. To put it concretely, we show how the num- 
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Figure 1: A pattern with a phyllotactic fraction (PF) a = 2/5. The 
angle between consecutive leaves (points) is 2/5 of a rotation. As 
a guide to the eye, a spiral, called the ontogenetic or genetic spiral, 
is drawn through the points in their order of succession. Repulsive 
interaction between the points is represented by (imaginary) prolate 
spheroids around the points. 
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Figure 2: Interaction V n (8 m —8 n ,m — n) between the points (m, 8 m ) 
and (n, 8 n ) is periodic in 8 m — 8 n with a period of 27r. By definition 
of X and n c , V n (8 m — 8 n ,m — n) vanishes for X < \8 m — 8 n \ /2tt < 1 
and \m — n\ > n c . 



ber of PTs, Apt, depends on the preset angle «o and the 
growth index n c (Fig. [5]). To the author's knowledge, no 
such mathematical evolutionary mechanism of phyllotaxis 
has ever been put forward. 

2. Model 

We restrict ourselves to the most common case of a 
spiral or helical pattern with a single 'leaf at each level 
(node). We use an integer n to label successive points 
along a genetic spiral (Fig. [1]), whose positions are given 
by the coordinates (9 n ,n) or by (x ni n) in terms of 9 n = 
2irx n . In this dimensionless representation with normal- 
ized length scales, the coordinate may be regarded to rep- 
resent the position of either leaf primordia on apical meris- 
tem or leaf traces in vascular system. Here the important 



point is that physical quantities are periodic with respect 
to the angular coordinate 9 n . The angle 9 n is measured 
from #o = for n = 0, and is regarded to take a value 
within — 7r < 9 n < n (or — 1/2 < x n < 1/2). To describe 
phenomcnologically a torsional force between two points 
(9 n ,n) and (9 m ,m) (n =/= to), we introduce repulsive in- 
teraction V n (9 m — 9 n ,m — n) (Fig. [2]). As a theoretical 
and phenomenological tool to clarify number theoretical 
structure of the macroscopic phenomenon, details for im- 
plementation of the interaction need not be specified. 
For the sake of simplicity and convenience, let us write 

V n {2ixx,m) = u n V{2irx,m) = u n v m V(x). (1) 

The angular dependence is represented by V{x) and the 
vertical (internode) dependence is described with v m , both 
of which are defined by the last equation. The factor u n 
represents the dependence on the subscript n of V n (2irx, m), 
which occurs because translational invariance along the 
stem (in the vertical direction) is broken in general. How- 
ever, it turns out that the factor u n drops out of our prob- 
lem. 

The interaction V n (9,m) is characterized by two pa- 
rameters. They are finite ranges of V n (9, to), in the angular 
direction 9 (or x — 9/(2tt)), and in the vertical direction 
to. For the former, we introduce a half-width X (< 1/2) 
of V(x), i.e., 

V{x) > 0, < \x\ < X, 
V{x) ~ 0, X < \x\ < 1/2. 

For example, let us use 

V(x)=e~(*)\ \x\<l/2. (2) 

Note that V(x) is periodic so that V(x) — V(x + 1), and 
we may set V^(O) = 1 arbitrarily. In terms of the width X 
thus defined, the angular width of the original interaction 
V n (9, to) is given by 

A0 ~ AttX. (3) 

The second parameter is the vertical range of influence n c . 
By definition, we have V n (9,m) > for < to < n c and 
V n (9,m) ~ for to > n c , or 

v n > 0, < n < n c , 
v n ~ 0, n c < n. 

In what follows, n c plays an important role. 

We investigate a regular spiral arrangement with the 
divergence angle 27ra, namely, 

9 n = 2irna. (x n — na.) (4) 

In principle, the phyllotactic index a may take any real 
number (Fig. [I}. Nevertheless, in our model, and in real 
life, a turns out to be a fraction, called phyllotactic frac- 
tion (PF). Phenomenologically, a obeys a mechanical re- 
laxation equation, 

da dE 
~di ~ ~~da' 



The right-hand side is a torsional force, represented by the 
total interaction E, 

oo 

E = Yl V ^ dm - n,m-n) (6) 

n— m>n 

oo oo 

= X! y »( 27rTOa > m )' ( ? ) 

n— m— 1 

Substituting Eq. ([I]), we obtain 

oo 

E = v(a)Yu n , (8) 

where 

oo 

v(a) = ^ VmV^rna). (9) 

The factor X^Ln Un m © ma y ^ e dropped hereafter, 
as it is a constant independent of a. Substituting Eq. (|5|) 
in Eq. ((SJ, we find the phyllotactic index a in static equi- 
librium as a minimum of the effective interaction v(a). 
For instance, when v(a) has a single minimum at a, we 
may assume a parabolic potential -E = | (a — a) 2 around 
the minimum a ~ a. Then we obtain ^ = —k(a — a) 
from Eq. ([5]), and a(t) = a + (a(0) — a)e as a solution. 
Therefore, in the end, we reach static equilibrium at the 
minimum a(oo) = a, irrespective of the initial value a(0). 
In general, v(a) may have many local minima. To which 
minimum a evolves into depends on the initial value a(0). 
This is the third parameter ao, 

a(0) = a . (10) 

We introduced several quantities to define our model. 
Among others, n c and ao are the most important. Main 
results given below are not affected essentially by the other 
quantities as X and v m for v(a) in Eq. (0). In the next sec- 
tion, we investigate conditions to realize a local minimum 
of v(a). Hereafter we restrict ourselves to < a < 1/2, 
because v(l — a) — v(— a) — v(a) by bilateral symmetry. 

3. Results 

3.1. Phyllotactic Fraction (PF) 

By way of illustration, v{a) in Eq. (j9]) is plotted for 
n c = 5,6 and 7 in Fig. [3j where we use Eq. ([2]) with X = 
0.1 and v n = 0.8™ for n < n c . 

Consider the case n c — 5, the solid curve in Fig. [3] 
We observe that v(a) has five minima. In effect, they lie 
around fractions a — 1/6,2/9,2/7,3/8 and 3/7. Let us 
increase n c from 5 to 7 to see if the minima are affected. 
For n c = 7, the two minima at 2/9 and 3/8 remain almost 
intact, whereas the other three minima 1/6, 2/7 and 3/7 
for n c = 5 are lost (Fig. |3J). 

To reach the minimum at PF 3/8, the initial value ao 
has to be in a range 1/3 < a < 2/5. Let us introduce 
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Figure 3: Effective interaction v(a) for n c = 5, 6 and 7 (X = 0.1 and 
v n = 0.8 n for n < n c ). To reach a minimum at a = 3/8, the initial 
value ao must fall within Act denoted by the double-headed arrow. 



the total width Aa of a range allowed for ao- For PF 3/8, 
we have Aa = | — 3 ~ 0.067. The range is indicated 
with the double headed arrow in Fig. [31 Similarly, we get 
Aa ~ 0.05 for 2/9, which is narrower than Aa ~ 0.067 for 
3/8. Therefore, if ao is to be chosen randomly for a fixed 
n c , it would be (0.067/0.05=1.3 times) easier to realize PF 
3/8 than to realize 2/9. 

From Fig. [3j we note that the minimum at 2/7 for 
n c — 5 and 6 turns into a local maximum for n c = 7. 
Similarly, the minimum at a = 3/8 becomes a maximum 
for n c — 8 (not shown in Fig. [3]), whereas the minimum at 
2/9 remains intact for n c — 5, 6, 7 and 8. For each PF, we 
define no and An by the condition 

no < n c < no 4- An (11) 

for the fraction to be a minimum. We obtain (no, An) = 
(5,2) for 3/8, and (n , An) = (5,3) for 2/9. 

When a minimum becomes a maximum, two new min- 
ima appear on both sides of the maximum. For example, 
the local minimum at a — 2/7 for n c — 5,6 ((no, An) = 
(4, 2)) becomes a local maximum at n c = 7. When a 
'mother' PF 2/7 is lost to become the maximum, it is 
flanked on both sides by two newborn 'daughter' minima 
at 3/11 and 3/10 (Fig. [3]). In other words, the mother PF 
a = 2/7 branches (fissions) into the daughter PFs 3/11 
and 3/10 at n c = 7. 

In this manner, we can locate the positions of the local 
minima of v(a) and construct their branching structure. 
They are shown in Fig. m where each dot representing a 
minimum is labeled with a fraction (PF). In Fig. 01 there is 
a vertical segment stretching upward from each dot. The 
length of the segment is An for the fraction. In Fig. [5l 
we plot the fractions in the An-Aa plane, as well as in 
the An-Aa-no space. In the figures, the observed main 
sequence of phyllotaxis is drawn with the bold line. 

In the above and the following, we make good use of 
the interesting and important property of the model that 
the divergence angle a of a phyllotactic pattern is given by 
a fraction (PF). Indeed, deviation from an exact fraction 
may be expected in general. This is discussed in Sec. 01 
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Figure 4: Position a of a local minimum of v(a) is represented by a 
point labeled with a fraction. Lines connecting the points represent 
hierarchical branching structure of the fractions. The length of the 
vertical segment stretching from a point represents An for the frac- 
tion (e.g., An = for 1/4 and An = 1 for 2/5). Fractions on the 
bold line comprise the main sequence of phyllotaxis. 
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a=2/ 5 (X=3) 
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Figure 6: An initial pattern (an = 0.382) of leaf primordia on the 
apex (top). A 'developing' leaf pattern for a = 2/5 (n c = 3) (middle) 
lies on top of a 'mature' leaf pattern for a = 1/3 (n c = 2) (bottom). 
Torsional motion is indicated by an arrow. 
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Figure 5: For lower order fractions in Fig. [4] Act is plotted against 
An. The inset shows a 3D plot of no against An and Aa (defined 
in Sec. 13. id . The bold arrows represent the main sequence. 



3.2. Phyllotactic Transition (PT) 

To compare with a real plant, we regard n c as a growth 
index. This is reasonable because the increase of n c in our 
normalized (cylinder) representation corresponds to a de- 
crease of plastochron ratio in a disc representation (fl6j|). 
and to a decrease of the internode distance if it were intro- 
duced explicitly ( [l[ ) . The index n c , presumably related to 
the length of leaf-trace primordia, is supposed to increase 
in the course of plant growth so that a stem is vertically 
organized into zones with different values of n c . In a devel- 
oping leaf zone near the apex, n c is large. In a mature leaf 
zone near the plant base, n c is small. For a given value of 
ao, we obtain patterns for different n c standing in a row 
along the stem, as illustrated in Fig. [5] Thus we explain 
phyllotactic transition (PT), the transition of a between 
different PFs along the stem. 

For given n c , from Fig. [5] we can read a resulting from 
arbitrary ao- Firstly, locate the reference point (aa,n c ) 
in the figure. Then, the PF a is selected from the neigh- 
boring two minima on both sides of the reference point 
by comparing the horizontal coordinate ao with that of a 
maximum between the two minima. If ao is on the left 
(right) of the maximum, then the minimum in the left 
(right) hand side is chosen. For example, for ao — 0.3, we 
get a = 1/2 for n c = 1, and a = 1/3 for n c = 2. For 
n c = 3, we choose a = 1/4 from the two minima 1/4 and 
2/5 because ao — 0.3 < 1/3, the maximum between the 
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Figure 7: Diagram for phyllotactic fractions (PF). PF a = 2/5 is 
obtained for n c = 3,4 and 1/3 < o < 1/2 (An = 1, Aa = 0.17). 
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two minima. 

As a result, we obtain a diagram in Fig. [7] The dia- 
gram may serve to get PF a for arbitrary pair (ao, rc. c ). In 
Table [TJ sequences of PF a for randomly chosen values of 
ao are given representatively. With real plants in mind, 
four of them are displayed in Fig. [8] 

3. 3. Mechanism of Phyllotaxis 

Nature's preference of particular PFs is accounted for 
simply by the following hypothesis: (H) A phyllotactic pat- 
tern of a fraction with larger An is more favored. 

For n c = 3 in Fig. 01 we have two possible patterns 
a = 1/4 or 2/5, depending on ao- According to (H), the 
latter 2/5 is favored, because An = 1 for 2/5 is larger 
than An = for 1/4. In this manner, as n c is increased, 
the hypothesis (H) gives us a sequence on the thick line in 
Figs. Hand El that is, 

1/2, 1/3,2/5, 3/8,5/13,- ■• . (12) 

This is the main sequence of phyllotaxis almost always 
observed. In fact, the main sequence covers more than 
90% of all observed cases ((HI). The numerator and the 
denominator of the PF in (TT2"j) are alternate numbers of a 
Fibonacci sequence, 

1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, 144, ■ • • , (13) 

in which every number is the sum of the preceding two. 
As a limit number of the sequence in (|12p . we obtain a = 
(3 - \/5)/2 ~ 0.381966 (27ra ~ 137.5°), the golden angle 
or the Fibonacci angle. 

We argue that the hypothesis (H) is biologically plau- 
sible, because the larger An ensures the more stability 
against expected variations of the growth index n c . To 
show this explicitly, let us define Npt as the number of 
PTs encountered in a sequence for a given ao. Npt counts 



Figure 8: Lower part of four sequences in Table are displayed 
vertically, (a) a = 0.146317 {N PT = 6), (b) a = 0.281978 (N PT = 
4), (c) a = 0.375791 (N PT = 3), (d) a = 0.469168 (N PT = 4). The 
most favorable is (c) with the least number of phyllotactic transitions 
(TVpx = 3). 



how many times PF a changes as the index n c is increased, 
so that it depends on ao and the upper bound of n c . Hence 
iVpT may be used as a measure of stability of a pattern 
with a given phyllotactic sequence. By definition, Npt 
is small for a sequence comprised of PFs with small An. 
Therefore, (H) may be rephrased as follows: (IT) A favor- 
able sequence has small Npt- 

Biological implications of (FT) may be understood in- 
tuitively from Fig. [5] According to (H'), the case (c) 
for ao — 0.38 (27rao = 137.5°) is the most favorable, 
because Npt = 3 for (c) is the smallest of all. What 
this means must be quite obvious from the figure. We 
give iVpx in the last column of Tabic [T] See the eighth 
row for a = 0.375791 in Table [TJ We find that A PT is 
the smallest for the main sequence, (fT2")l . for ao ~ 0.38 
(27rao = 137.5°). As a second sequence with small Npt, 
we notice a sequence for ao — 0.28 (27ra = 99.5°), the 
fifth row in Table [1] and (b) in Fig. [5] This sequence is 
observed but less commonly, and sometimes called the first 
accessory sequence. In general, the main sequence always 
sets the lower limit of Npt, although there may be other 
sequences with the same lowest value. 

In Fig. O iVpx for n c up to 4, 5, 6, and 7 is plotted 
against ao- The figure indicates how the samples with 
different ao are discriminated as they grow. In the pro- 
cess of increasing n c to 4, any sample for 1/3 < a < 1/2 
(Aa = 0.17) is more favorable than that for < ao < 1/3, 
because the former has a smaller value JVpx = 2 than the 
latter with Npt = 3. As we increase n c to 7 (the solid line 
in Fig. |9]), only restricted samples within a narrow win- 
dow 1/3 < a < 2/5 (Aa = 0.07) are favored because of 
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Table 1: Sequences of PF a resulting from arbitrarily chosen initial values of ao. TVpt in the last column counts the number of phyllotactic 
transitions (the number of times PF a changes) for 1 < n c < 17. 
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Figure 9: The number of phyllotactic transitions iVpx for n c = 4, 5, 
6, and 7 are drawn against ao- As n c increases to 7, initial values ao 
within 1/3 < ao < 2/5 (Aa = 0.07) are specifically selected because 
iVpT is the smallest there. 



iVpx = 3, which is the lowest value for n c — 7. In this fil- 
tering process, the width Aa for ao decreases rapidly as n c 
increases. In Fig. [101 we plot Npt against ao and n c . The 
figure indicates that the minimum of Npt with a decreas- 
ing width Aa develops around ao — 0.38 (27rao = 137.5°) 
as n c increases. In the end, a single value is selected for 
ao, namely, the golden angle ao = (3 — v / 5)/2 ~ 0.381966 
(Aa ->■ 0). 

In Fig. 1111 we plot iVpT for n c up to 20 and 50, indicat- 
ing how the stepwise ridges and troughs of Apt develop 
as n c increases. As n c increases, ao should be gradually 
trapped in a minimum of Ap-p- We find a small and shal- 
low local minimum around ao — 0.22 (27rao = 78°), be- 
sides the absolute minimum at ao ~ 0.38 and the first ac- 
cessory ao — 0.28 mentioned above. These three cases are 
collectively called normal phyllotaxis (Sec. l3.5l and (B.13)). 
It is remarked that Apt (or Apt / n c ) becomes singular in 
the limit n c — > oo, for all the steps are subdivided and the 
widths of the steps shrink without limit as n c increases. 
This is clear from Fig. rT2J in which n(ao) = Npt/ti c is 




Figure 10: JV PT as a function of ao for 1 < n c < 20 (cf. Fig. [9J . A 
narrow trough develops around ao ~ 0.38 (27rao = 137.5°), as n c 
increases. 



plotted for n c =100 and 500. Note that n(ao) is mini- 
mized at the golden angle. But, from the figure, it appears 
almost impossible to reach the golden angle continuously 
by variational optimization. In our static mechanism, the 
golden angle is singled out through the screening process. 

In Figs. Q2] and [TU the range of ao within which the 
number Apt is the smallest for given n c is filled with the 
solid horizontal bar. It is clearly shown that several values 
are specifically favored for ao. Among others, we remark 
that the golden angle 27rao — 137.5° (a ~ 0.382) is sin- 
gled out for n c near but less than the Fibonacci numbers 
in (TI"3")) . The accuracy Aa for the golden angle converges 
rapidly, as mentioned above. We have Aa = 1/15 for 
n c = 7 in Fig. |3 and Aa = 1/1870 around n c ~ 85 (cf. 
Table lIOll. 

To conclude, we may regard the hypothesis (H) as a law 
of phyllotaxis. In plain words, the Fibonacci phyllotaxis 
with the main sequence is favorably singled out because of 
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Figure 11: iVpT for n c = 20 (solid line) and 50 (dashed line). The 
arrow at the minimum of Np^ indicates the golden angle ao ~ 0.3820 
(2na = 137.5°). 





Figure 13: The range of ao within which Wpx takes the smallest 
value for given n c is filled with horizontal bars, thereby favorable 
values for ao are pointed by 'tapering needles'. They are ao ~ 
0.38,0.28,0.22 (normal phyllotaxis), a ~ 0.42,0.44, and also a ^ 
0.30, 0.37 (cf. Sec.E}. The golden angle a ~ 0.38 ( 27ra = 137.5°) 
is singled out for n c ~ 20, 30, 50, 80. 



Figure 12: NpT/n c for n c = 100 (dashed line) and 500 (solid line). 
The minimum at ao ~ 0.38197 is indiscernible for n c = 500 (Ao ~ 
10~ 5 ). 



its special stability against inevitable structural changes 
expected in the growing process. 

3.4- Mathematics 

The main results presented above are obtained directly 
from simple numerical analysis as shown in Fig. [3l Indeed 
it is straightforward to check them by hand when n c is 
small, but it soon gets complicated as n c becomes large. 
Mathematical analysis helps us not only to derive useful 
formula but also to deepen understanding of the mathe- 
matical structure of the problem. In particular, it is help- 
ful for us to have analytical expressions of (no , An) and 
Aa for PFs belonging to typical sequences found in na- 
ture. The mathematical analysis is indispensable if we do 
not content ourselves with several specific circumstances of 
lower phyllotaxis. Mathematical results are delegated to 
the appendixes. Here we point out only that Fig.[4]has the 
same structure as the Stern-Brocot tree of number theory 
which contains each rational numbers exactly once. 
Relations between numbers in the tree are concisely repre- 



100 




0.39 



Figure 14: The range of ao for the least A^px is filled with horizontal 
lines (Fig. [131 enlarged vertically and horizontally). 
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sented in terms of mediants ( |Appcndix A I and continued 
fractions ( |Appendix B 1 . 



3.5. Sequences 

As given in Table [TJ a general phyllotactic pattern 
derived from a random value of ao does not fit in with 
observed regularity. Indeed, if PT would have to occur 
too frequently, the concepts of the PF and PT themselves 
could become indefinite. As a matter of fact, fortunately, 
observed sequences come around with astonishing regu- 
larity. Only several types of sequences exist in nature. 
For the purpose of classifying the observed sequences, one 
often adopts a tacit theoretical procedure of inferring a 
mathematical limit a of a given sequence by extrapolation. 
The limit divergence a, generally an irrational number, is 
then used to represent the sequence. However, it must be 
kept in mind that any phyllotactic sequence does terminate 
finitely in practice, and inferring a limit from a finite se- 
quence may be problematic. Be that as it may, a sequence 
of principal convergents of a noble number (Eq. (|B.9p ) has 
been a central subject. We can assess frequencies of oc- 
currence of various sequences quantitatively by compar- 
ing iVpT for the sequences. Below we use a shorthand 
bracket notation for a sequence given in a paragraph be- 
low Eq. (|B.14p in |Appendix B 



The main sequence with the limit divergence ao — 
0.382 (27ra = 137.5°) is given by 

[2] : 1/2, 1/3, 2/5, 3/8, 5/13, 8/21, 13/34, 21/55, • • • . (14) 

The first accessory sequence with the limit ao = 0.276 
(27ra = 99.5°) is 

[3] : 1/2,1/3, 1/4,2/7,3/11,5/18,8/29, 13/47,- •• . 

The second accessory sequence for ao = 0.217 (27rao = 
78.0°) is 

[4] : 1/2, 1/3, 1/4, 1/5, 2/9, 3/14, 5/23, 8/37, 

These are called normal phyllotaxis. As an example of 
anomalous phyllotaxis, the first lateral sequence for ao = 
0.420 (27ra = 151.1°) is 

[2, 2] : 1/2, 1/3, 2/5, 3/7, 5/12, 8/19, 13/31, 21/50, 

In addition, we can think of the sequence for ao = 0.296 
(27ra = 106.4°), 

[3, 2] : 1/2, 1/3, 1/4, 2/7, 3/10, 5/17, 8/27, 13/44, 

And, for a = 0.367 (27ra = 132.2°), 

[2,1,2] : 1/2,1/3,2/5,3/8,4/11,7/19,11/30,18/49, ••• . 

There is controversy concerning the existence of the last 
two sequences ( [27 , 11 [ ) . As the general fact of observation, 
any other sequence than the main sequence [2] may be 
regarded as exceptional. 




Figure 15: TVpx is plotted against n c for major sequences. The 
sequences with the limit index a ~ 0.38,0.28,0.22,0.42,0.30 and 
0.37 are represented as [2], [3], [4], [2,2], [3,2] and [2,1,2], respectively. 
The main sequence [2] sets the lower limit of iVpx- 



For these sequences, Apt is plotted against n c in Fig.[T5l 
The main sequence always sets the lower limit of A^px- At 
n c = 4, [2,2] and [2,1,2] have Af PT = 2, and the others 
have A^pt = 3. This is because the first three terms of 
[2,2] and [2,1,2] (up to 2/5) are the same as the main se- 
quence [2]. From around n c = 10, the priority order of 
occurrence of the sequences is inferred as [2,2], [3], [3,2], 
and [4]. For such small n c as shown in Fig. 1151 the order 
of [2,1,2] and [2,2] cannot be decided uniquely. Note that 
[3,2] and [2,1,2] may be regarded as satellite sequences of 
[3] and [2] , respectively. This is seen from Fig. [T3] There- 
fore, according to our result, extraordinary [3,2] and [2,1,2] 
are less unlikely than [4] of 'normal' phyllotaxis. In Ta- 
ble [21 for general sequences, we tabulate the values of n c 
at which A^pt is minimized. Many sequences become fa- 
vorable when n c is equal to a Fibonacci number in (|13[) . 
This is clear from Fig. [T31 We find that normal phyllotaxis 
of higher order [n] (n — 4, 5, ■ • • ) appears not so specially 
preferable as widely supposed. In fact, Fig. H3l indicates 
that Afpx has no absolute minimum for < ao < 0.2 and 
n c > 4 (a = 0.178 for [5] and a = 0.151 for [6]). A 
low-order sequence with [l,m,n] (l,m,n < 3) looks rather 
favorable. 

The order in the frequency of occurrence is generally 
consistent with the available data, though the number of 
observations is still too limited to draw a definite conclu- 
sion (@,[H|). 

4. Discussion 

We have made full use of an important result of our 
model that the phyllotactic index a is given by a fraction. 
As a problem of macroscopic physics, it goes without say- 
ing that this is but a good approximation and mathemat- 
ical rigor should not be expected in this respect. In this 
section, we investigate the effect of the lateral width X of 
the interaction V(x) in Eq. ([2]). There is an optimal value 
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Table 2: Table for the growth index n c at which a given sequence makes a pattern with the minimum TVpx- In the first row, the minimum 
value TVpx is given. The main sequence [2]: (1,2) covers all integers for n c (iVpx = 2 for n c = 3,4, etc.). The bracket notation for the 
sequence in the first column is explained in a paragraph below Eq. IIB.14II . (A pair of numbers in the round bracket are parastichy pairs, i.e., 
seed values for a Fibonacci recurrence relation. ) 



Npt 


2 


3 


4 


5 


6 


7 


8 


9 


[2] : 


(1,2) 


3,4 


5-7 


8-12 


13-20 


21-33 


34-54 


55-88 


89-143 


[3]: 


(1,3) 


3 


5,6 


8-10 


13-17 


21-28 


34-46 


55-75 


89-122 


[4] : 


(1,4) 


3 




8 


13 


21,22 


34-36 


55-59 


89-96 


[5] : 


(1,5) 


3 
















[2,2' 




(2,5) 


3,4 


5,6 


8-11 


13-18 


21-30 


34-49 


55-80 


89-130 


[2,3 : 




(2,7) 


3,4 


5,6 


8 


13-15 


21-24 


34-40 


55-65 


89-106 


[2,4 




(2,9) 


3,4 


5,6 


8 












[3,2[ 




(3,7) 


3 


5,6 


8,9 


13-16 


21-26 


34-43 


55-70 


89-114 


[3, 3 




(3,10) 


3 


5,6 


8,9 




21,22 


34,35 


55-58 


89-94 


[3,< 




(3,13) 


3 


5,6 


8,9 












[4,2" 




(4,9) 


3 




8 




21 


34 


55,56 


89-91 


[4,3 




(4,13) 


3 




8 












[2,1, 


2] 


(3,8) 


3,4 


5-7 


8-10 


13-18 


21-29 


34-48 


55-78 


89-127 


[2,1, 


3] 


(3,11) 


3,4 


5-7 


8-10 


13 


21-24 


34-38 


55-63 


89-102 


[2,1, 


4] 


(3,14) 


3,4 


5-7 


8-10 


13 










[2,2, 


2] 


(5,12) 


3,4 


5,6 


8-11 


13-16 


21-28 


34-45 


55-74 


89-120 


[2,2, 


3] 


(5,17) 


3,4 


5,6 


8-11 


13-16 


21 


34-38 


55-60 


89-99 


[2,2, 


4] 


(5,22) 


3,4 


5,6 


8-11 


13-16 


21 








[2,3, 


2] 


(7,16) 


3,4 


5,6 


8 


13-15 


21,22 


34-38 


55-61 


89-100 


[2,3, 


3] 


(7,23) 


3,4 


5,6 


8 


13-15 


21,22 








[2,3, 


4] 


(7,30) 


3,4 


5,6 


8 


13-15 


21,22 








[2,4, 


2] 


(9,20) 


3,4 


5,6 


8 












[3,1, 


2] 


(4,11) 


3 


5,6 


8-10 


13,14 


21-25 


34-40 


55-66 


89-107 


[3,1, 


3] 


(4,15) 


3 


5,6 


8-10 


13,14 










[3,2, 


2] 


(7,17) 


3 


5,6 


8,9 


13-16 


21-23 


34-40 


55-64 


89-105 


[3,2, 


3] 


(7,24) 


3 


5,6 


8,9 


13-16 


21-23 








[3,3, 


2] 


(10,23) 


3 


5,6 


8,9 




21,22 




55 




[3,3, 


3] 


(10,33) 


3 


5,6 


8,9 




21,22 








[4,1, 


2] 


(5,14) 


3 




8 


13 










[4,1, 


3] 


(5,19) 


3 




8 


13 










[4,2, 


2] 


(9,22) 


3 




8 




21 









X for X around which a is given by a fraction, as expected. 
The optimal width X may be used to see if a minimum of 
v(a) is really achieved practically. If the potential function 
v (a) is nearly constant and flat in a wide region around a 
minimum, it will take too long time to reach the true min- 
imum because the torsional driving force — dv j°^ toward 

da 

the minimum, in the right-hand side of Eq. ([5]), becomes 
practically zero. This happens if X <C X, as shown in 
Fig. [T5]for X = 0.05. It is known that the use of a frac- 
tion, as originally made by Schimper and Braun, is not 
always adequate ([!]). 

To put it concretely, let us examine PF 3/8, which 
occurs in between 1/3 and 2/5. We consider 

l/3<a<2/5, (15) 

and 

5<n c <8. (16) 



Then, for v(a) in Eq. ([9]), all the other terms than m = 3 
and to = 5 are effectively neglected, so that we get 

v(a) ~ v 3 V(3a) + v 5 V(5a) 

= v 3 V(3a~l)+v 5 V(5a-2), (17) 

for V{x) is periodic. In effect, this is a minimal model 
for PF 3/8. This may be interpreted as a mathematical 
expression of Hofmeister's rule, that new 'leaf arises in 
the largest gap between the previous ones. Note, however, 
that we are not concerned with mechanical dynamics of 
phyllotaxis at the apex. 

On the one hand, as a function of a, y(3a — 1) has a 
peak with the width X/3 at a = 1/3, the lower boundary 
of ([T5|) . On the other hand, V(5a — 2) has a peak with 
the width Xjh at a = 2/5, the upper boundary of (fl~5|) . 
An optimal width X is estimated by equating the total 
width X/3 + X/5 with the allowed range for a, that is, 
Aa = 2/5 - 1/3. We obtain the result X = 1/8 for 3/8. 
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1/3 2/5 



a 

Figure 16: For various X around an optimal value X = 0.125, the 
effective potential v(a) is plotted by using Eq. J2J and v n = 0.6 n for 
n < n c (5 < n c < 8). For X = 0.05 (< X), v(a) has a flat bottom 
for a wide range of a. 



5. Conclusions 

With the aid of biological hypotheses, we showed that 
prevalent sequences of phyllotactic fractions are satisfac- 
torily explained by a physical model of plant growth. The 
model has interesting mathematical properties. Among 
others, we bring to light a Stern-Brocot type number- 
theoretical structure that has been unnoticed thus far in 
this field. To extract the mathematical essence of the 
phenomenon, we have to base our theory on the abstract 
model by discarding real biological implementation as non- 
essential details. According to the proposed static mech- 
anism, the phyllotactic pattern with the main Fibonacci 
sequence is naturally selected because it entails the least 
number of phyllotactic structural transitions while grow- 
ing to a mature plant. 



According to Eq. (J3J , the optimal angular width of inter- 
action is A9 ~ 90° for 3/8. In the ideal case X — X, v(a) 
has a minimum ata = l/3 + A/3 = 2/5 — A/5 = 3/8, as 
expected. The derivation outlined here indicates that the 
result will not depend on V(a) specifically. 

We present Fig. \TE\to show the effect of X on v(a). To 
draw the figure, we use a full form of v(a) and did not use 
the approximation in Eq. (|17[) . In any case, the difference 
due to the approximation is negligible. As mentioned, for 
X ~ X, we obtain a ~ 3/8 properly with good accuracy. 
As X is decreased from X, v{a) around the minimum gets 
flattened. As a matter of practical fact, the minimum 
would not be reached in the extreme case of A <C X. 
When X < X, there appears a flat region in 



IX 2 X 

3 + J <a< 5~J> 



(18) 



the width of which is Aa flat = (1 - 8A)/15 = (1 - f )/15. 
This is smaller than the full width Aa = 1/15 of (fTS]) , as 
it should be. In practice, insofar as a stays within the flat 
region, we may rather observe the primary angle a — ao 
as it is, since the secondary torsion is not effective any 
longer. 

In general, the optimal width X for a fraction p/q is 
given by X — 1/q, and the width of the flat region is 



An lbli = ( 1 - J ) An " 



(19) 



for X < X ( | Appendix C[ ). As we follow any sequence 
up in the tree of Fig. HJ the denominator q n of PF a = 
Pn/On stays constant or increases, so that X is constant or 
decreases. In effect, it decreases roughly as X ~ l/n c . By 
contrast, the model parameter X should be fit with a real 
plant. Consequently, the condition X < X may hold in 
early few terms of a sequence (namely, 1/2, 1/3, etc., when 
n c is small). Then we could no more expect to observe 
these low order PFs. We rather observe an inherent value 
a = ao so far as it falls within a flat region with the width 
Aaflat. This is consistent with observations. 
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Appendix A. Mediant 

The mediant of two fractions m/n and p/q is given by 
(m + p)/(n + q), where m and n (p and q) are relatively 
prime integers. Let us call m/n and p/q as parent fractions 
of a child fraction (m + p)/ (n + q). 

The Stern-Brocot tree of fractions between and 1/2 is 
obtained by the following operation (Q). Start from the 
initial fractions 0/l(=0) and 1/2. Repeat inserting the 
mediant of two adjacent fractions, and arranging them in 
numerical order. The first mediant is 1/3 between 0/1 
and 1/2, and they are arranged as (0/1), 1/3, (1/2). In 
the second order, we obtain 1/4 between 0/1 and 1/3, and 
2/5 between 1/3 and 1/2. They are arranged as (0/1), 
1/4, (1/3), 2/5, (1/2). In the third order, we obtain (0/1), 
1/5, (1/4), 2/7, (1/3), 3/8, (2/5), 3/7, (1/2). Here we 
put the fractions in the previous orders in parentheses. 
The Stern-Brocot tree has been anticipated by Schimper 

(0,113). 

Our tree in Fig. [4] is related to but not the same as the 
Stern-Brocot tree. As an important difference, we have to 
order fractions by the growth index n c . In other words, we 
need to know (no, An) of the fractions. For instance, let us 
consider 2/5. For n c = 3 in Fig. HI (m+p)/(n + q) = 2/5 
is the mediant oip/q = 1/3 and m/n = 1/2 (m = l,n = 
2,p = 1, q = 3). By ordering the fractions according to the 
denominators, q and n, let us call p/q = 1/3 and m/n = 
1/2 as the younger and the older parent of 2/5. On the one 
hand, the mediant 2/5 is born when the younger parent 
p/q = 1/3 dies at n c = q = 3. On the other hand, 2/5 
dies at n c = n + q = 5, the denominator of 2/5. Therefore, 
we obtain 3 < n c < 5 and An = 1 for 2/5. In general, 
the mediant (m + p)/(n + q) of m/n and p/q (n < q) has 




0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

a 

Figure A. 17: The branching structure in Fig. [3] is schematically 
shown (a) for p/q < m/n and (b) for p/q > m/n, where m and n, p 
and q are relatively prime numbers satisfying n < q. In both cases, 
two daughter fractions (m+2p)/(n + 2q) and (2m+p) / (2n + q) at the 
top of the figure are derived from a mother fraction (ni + p)/ (n + q) 
in the middle. The mother (m+p) /(n + q) is derived as a child of the 
parents p/q and m/n. We obtain An = q — 1 for (m + 2p)/(n + 2q), 
and An = n — 1 for (2m + p)/ (2n + q) and (m + p) / (n + q) . The 
path to increase An traces the bold zigzag line. 



q < n c < n + q. For (m + p)/(n + q) (n < q) , we obtain 

(n ,An) = (q,n-l). (A.l) 

Next we consider Aa. In order to realize a = (m + 
p)/(n + q), ao has to lie between the parents m/n and 
p/q, and we get Aa = |~ — H|. The parent fractions 
belonging to the tree arc shown to satisfy \mq — np\ = 1 
([§]). Hence we conclude 

Aa=— (A.2) 
nq 

for (m+p)/ (n + q). These formulas may be used for Fig. [5) 
For the hypothesis (H) in Sec. 13.31 we have to compare 
An of two daughter fractions derived from a mother frac- 
tion in Fig. SI (We tell a mother-daughter relation from 
a parent-child relation.) Consider (m + p)/(n + q) as a 
mother fraction, derived as a child of parents m/n and 
p/q (n < q) (Fig. IA.17)) . One daughter fraction (2m + 
p) / (2n + q) occurs between (m+p) /(n + q) and m/n. The 
other daughter fraction (m + 2p)/(n + 2q) occurs between 
(m+p) /(n + q) and p/q. On account of Eq. (|A. 1|) . the for- 
mer (2m+p)/(2n+q) has (no, An) = (n+q, n— 1), whereas 
the latter (m+2p) / (n+2q) has (n ,An) = (n+q, q— 1). By 
assumption n < q, the daughter fraction (m + 2p) / (n + 2q) 
has the larger An = 17—1 than (2m + p)/ (2n + q) with 
An = n — 1. For two daughter fractions derived from a 
mother fraction, the one with a larger denominator always 
increases An, whereas the other does not change An from 
the mother fraction. 

This rule is read from Fig. [5] At every branching point 
(node), one branch grows to the right to increase An. The 
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io\ 




5.9) 



(7,9 1 



(4,7) 



(3,7) 



(3,8) 



(5.8) 



-(6,7) (5,6) 



(1,7) 



(4,5) 



- (5,7) (2,7) 



(2.5I 



Figure A. 18: Patterns with the same divergence angle a = 0.38 
appear to have different parastichy pairs (n, q) = (3,5) (left) and 
(5,8) (right), depending on a vertical length scale. The 5-parastichy 
winds up clockwise. Thus this 'transition' of a parastichy pair is a 
superficial phenomenon. The apparent change from (3,5) to (5,8) is 
called rising phyllotaxis. Rising phyllotaxis restricts a range allowed 
for a (Fig. [09t . 



(2.3: 



other goes down along the ordinate in the main figure. Ac- 
cording to (H), the sequence with ever increasing An com- 
prises the most favorable branch of our evolutionary tree. 
A favored sequence in the tree diagram of Fig. 2] traces 
a zigzag path as depicted with the bold line in Fig IA. 171 
This result may be compared with a dynamical counter- 
part ( |Appcndix B| [Appendix Cj ). 

In practice, when the divergence angle a deviates from 
an exact fraction, a pair of integers (n,q), called a paras- 
tichy pair, is used to represent a phyllotactic pattern con- 
sisting of most conspicuously visible families of n and q 
secondary spirals (parastichies) crossing with each other. 
A g-parastichy is a secondary spiral running through the 
points (p, 9 P ), (p + q, p+q ), (p + 2q, 6> p+ 2g ), (p + 3q, 6 p+3q ), 
etc., where p = 0, 1, 2, • • • , q - 1 (Fig. IA.18[t . To derive a 
parastichy pair for given a is a purely geometrical prob- 
lem (HQ). In our tree system, we find a simple result 
that we obtain a visible parastichy pair (n, q) when a lies 
between two neighboring (parent) fractions m/n and p/q 
in our tree. Accordingly, PF a = (m + p)/(n + q) may be 
replaced by the parastichy pair (n,q). In place of Fig. [JJ 
we obtain Fig. IA. 191 which may be useful to analyze real 
systems. 



Appendix B. Continued Fraction 

A real number a (0 < a < 1) is represented as a con- 
tinued fraction, 



1 



a 



[ai,a 2 ,a 3 , 



(B.l) 



«i 



0-2 



«3 



where di (i = 1, 2, • • • ) is an positive integer. Every ratio- 
nal number has two continued fraction expansions. In one 
the final term is 1, that is, a rational number a is repre- 
sented finitely as a = [<Zi,a2,-- - , a n , 1], For an irrational 
number a, there is a successive rational approximation, 
Pn/qn = [fflij • • • j a-n] ( n — 1,2,---), in terms of relatively 



Figure A. 19: Diagram of a phyllotactic system with one genetic 
spiral, showing the visible opposed parastichy pair (n, q). For 
3/8 < a < 2/5, we may obtain (1,2), (2,3), (3,5) or (5,8), depending 
on a vertical length scale (Fig. I A. 1811 . Conversely, if we observe ris- 
ing phyllotaxis through parastichy pairs (1,2), (2,3), (3,5) and (5,8), 
then the divergence angle a must meet the constraint 3/8<a<2/5 
(Act = 0.025). For a system with J genetic spirals, a and (n, q) are 
replaced by Ja and J(n, q), respectively. 



prime positive integers p n and q n satisfying the recursion 
relations, 



and 



Pn+2 = a n+2 Pn+l+Pn, 
q-n+2 — CLn+2qn+\ + 9nj 

Po = 0, pi = 1, q = 1, qi = a\. 



(B.2) 
(B.3) 

(B.4) 



The fraction p n /q n (n = 1, 2, ■ ■ • ) is called a principal con- 
vergent of a. The difference between successive principal 
convergents is given by 



Pn+l 



Pn 



(-iy 



Qn+l In <Zn<7n+l 
and a lies between even and odd order convergents, 

P2k ^ ^ P2k+1 

— < a < . 



(B.5) 



q2k 



92fe+l 



(B.6) 



Thus the principal convergent p n /q n approaches to the 
limit a in a zigzag manner. 

In terms of the bracket notation defined in Eq. (IB. 1|) . 
we obtain Fig. IB.20l in place of Fig. 2) From the figure, we 
immediately notice the bifurcation rule holding at every 
node, 



[<2l, 02, • • • ,a n , 1] 



[ai,a 2 , 
\ai,a 2l ■ 



■ ,a r 



,1,1] 
-1,1]. 



(B.7) 



It is easily checked that the upper branch in (|B.7|) increases 
An, and the lower one conserves An d Appendix A I. There- 
fore, according to (H) in Sec. 13.31 it is particularly impor- 
tant for us to study the following sequence of fractions 
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[1.1] 



0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

a 

Figure B.20: Fractions in Fig. [4] are represented with the bracket 
notation for continued fractions defined in Eq. IB.lt . 



(<Zi ^ 1 for i = 1, 2, • • • , n): 

[ai,a 2 , ■■■ ,a n , 1], 

[ai,a 2 , • • • ,a„,l, 1], 

[ai, o 2 , • • • , a n , 1, 1, 1], 

[ai,a 2 , • • • ,a„, 1, 1, 1, 1], • • • . (B.8) 

These are comprised in the principal convergents of a noble 
number, 

"noble = [ai, 02, ' ' ' , On, 1, 1, 1, 1, " •]• ( B - 9 ) 

Owing to the succession of a, = 1 [i > n), for a sequence 
of fractions p n jq n on a favored branch according to (H), 
the numerator p n and denominator q n obey the Fibonacci 
recursion relations p n — p n —i + Pn—2 an d In = Qn—l + 
q n -2- 

A sequence of principal convergents of a noble num- 
ber has been given a special status in theoretical stud- 
ies of phyllotaxis (0, [H, [HI [l^]). The most important 
is the golden angle (per 2tt) for a\ = 2 and a n = 1 
(n = 1,2,3,---). 

a = [2, 1, 1, 1,. •-] = 1/(2 + r" 1 ) = t- 2 = (3 - V5)/2, 

(B.10) 

where the golden ratio r is given by 

t=1+ [1,1,1," •] = (>/5 + l)/2 ) (B.ll) 

or 

r- 1 = [l,l,l,».] = (>/5-l)/2. (B.12) 

In the literature, a sequence is referred to in several 
ways. In particular, a sequence with a limit index 

a=[i,l,l,l,..-] = l/(t + r- 1 ) (B.13) 



(t = 2,3, • • • ) is called normal phyllotaxis. The main se- 
quence for Eq. (jB.lOp corresponds to the special case t = 2. 
The cases t — 3 and t = 4 are called the first and the sec- 
ond accessory sequence, respectively. On the other side, a 
sequence with a limit angle 

a= [2,t, 1,1,1, •••] = 1/(2 + (t + r- 1 )- 1 ), (B.14) 

(t = 2, 3, • •• ) is sometimes called the lateral sequence. 

For convenience sake, let us introduce another concise 
notation. To denote a whole sequence with a limit of 

a = [ax,a 2 , ■ ■ ■ , a n , 1, 1, • • • ] = l/(ai+(a 2 +( h(a„+T _1 ) _1 ) _ 

(B.15) 

we simply reuse a symbol for a continued fraction [ai, a 2 , • • ■ , a„ 
(a n 7^ 1). Hence the main sequence is 

[2] : 1/2,1/3,2/5,3/8,5/13,8/21, 13/34,21/55,-" . 

(B.16) 

To obtain the sequence [3, 2] with a limit a = 1/(3 + (2 + 
I— 1 )" 1 ), collect the fractions from 1/2 to [3,2] = 2/7 
up along the tree in Fig. |4j namely, 1/2, 1/3, 1/4, 2/7. 
Then, from the daughter fractions 3/10 and 3/11 of 2/7, 
choose the unfavorable one with a smaller denominator, 
namely, 3/10. After them follow all the fractions on the 
favorable branch ramifying from 3/10, namely, 5/17, 8/27, 
13/44, 21/71, 34/115, 55/186, 89/301, • • • . As noted be- 
low Eq. (|B.9|) . these obey the Fibonacci recursion relation 
(e.g., 13/44 = (5+8)/(17+27)). To summarize, we obtain 

[3, 2] : 1/2, 1/3, 1/4, 2/7, 3/10, 5/17, 8/27, 13/44, 21/71, 34/115, 

Major sequences are given in Sec. 13.51 

A sequence is commonly represented by parastichy num- 
bers instead of fractions. Translation is made without dif- 
ficulty by noticing the denominators ( | Appendix A[ ) . 

[2] : (1, 1), (1, 2), (2, 3), (3, 5), (5, 8), (8, 13), (13, 21), 

(B.17) 

[3,2] : (1,1), (1,2), (1,3), (3, 4), (3, 7), (7,10), (10,17), (17,27),- 

(B.18) 

We may combine repeated numbers on favored branches. 

[2] : 1,2,3,5,8, 13, 21, • • • . (B.19) 

[3, 2]: (1,2), (1,3, 4), (3, 7, 10,17, 27,--. (B.20) 
All but the last branch may be omitted. 

[3,2] : 3,7,10,17,27,- •• . (B.21) 

Most concisely, only the first two numbers may be given 
as seed values of recurrence. 

[2]: (1,2). (B.22) 

[3,2]: (3,7). (B.23) 

In this notation, typical sequences are generally given as 
follows (cf. Table©. 

[a]: (l,o). (a = 2, 3,-".) (B.24) 



\p,a]:{p,ap+l). (a,p = 2, 3, • • • .) (B.25) 

\p-l,l,a-l]:(p,ap-l). (a,p = 3,4,--- .) (B.26) 

[p-l,2,o-l] : (2p-l,o(2p-l)-p). (o,p = 3,4,--- .) 

(B.27) 

According to 0], the first sequence system [a] is common, 
the second sequence system [p, a] is rarely observed, and 
the third sequence system \p— 1, 1, a— 1] is extremely rare. 
This is consistent with our results in Table [21 

Now our task is to get Aa and (no, An) for a se- 
quence of principal convergents of a noble number. From 
Eq. (|B~5|) . we get 



Pn 
In 



Pn-l 

q n -i 



(-iy 



qn-lQn 



and 



Pn_ _ Pn-2 _ (-l)"a„ 



(B.28) 



(B.29) 



?n-2 QnQn-2 

These equations signify that p n / q n lies between p n -i/ q n - 1 
and p n -2/q n -2- We find that p„/g n has 



Aa 



Pn-l Pn-2 



q-n-l q-n-2 



qn-iqn-2 



(B.30) 



According to our model, a child fraction p n /q n is born at 
n c — q n -i from a parent fraction Pn-i/q-n-i, and dies at 
n c = q„. Henc e, we obtain g„_i < n c < q n for p n /q n - 
Using Eq. (IB. 3ft . for p n /q n of a noble number with a n = 1, 
we conclude 



(n , An) = (<7„„i,(j„- 



1) 



(B.31) 



To put it more concretely, hereafter we restrict our- 
selves to the most important case of the golden angle 
in Eq. (|B.10[) . The principal convergent p n /qn, Aa and 
(no, An) are presented in Table lB~3l In this simplest case, 
we have p n = F n and q n — F n+ 2 , where F n is the Fibonacci 
number defined by the recurrence 



= 

Fi = 
F n = 



0, 
f, 

F n —i 



F n . 



(n>l) 



(Table lBT4l ) The number of phyllotactic transition A^pt 
simply counts the number of p n / q n . Consequently, we ob- 
tain iVpT = 2n + 1 for 



P2n ^ ^P2 

< a < 

<72n 



n+1 



92n+l 



(B.32) 



and 



<72n+i < n c < q 2n +2- (B.33) 

In terms of F n , we obtain Npt = In + 1 for F2 n +3 < n c < 
F2n+4, and Aa = I / '(i*2n+2^2n+3) • As a result, irrespec- 
tive of whether n is even or odd, we obtain 



Apt = n 



(B.34) 



for 
and 



F n +2 < n c < F„ + 3, 



Aa 



F n +iF n +2 



Using a formula ([s| 
f 



(B.35) 
(B.36) 

(B.37) 



we may regard n c ~ r A ' PT+2 - 5 /v / 5 by (|B~35|) . Then, we 
get 

iVpT = log(v / 5n c )/logr - 2.5, 



and 



Aa = (r/n c 



(B.38) 



(B.39) 



Thus Npt/tIc and Aa vanish in the limit n c — > oo. The 
logarithmic dependence in Eq. (IB.38[) for the irrational 
number a = t -2 is contrasted with a linear dependence 
Npt — n c /n for a rational number a — 1/n (Fig. I12j) . 
These results are used to guess a growth index n c inversely. 
To achieve accuracy of Aa = 0.02 (i.e., a = 0.38 ± 0.01), 
we need n c ~ 11 according to Eq. (|B.39|) . To obtain a 
parastichy pair (-FiO)-Fii) = (55,89) of a sunflower head, 



we have to attain Aa = 



1/4895, for which we 



FioFu 

need n c ~ 113. Indeed this lies between Fn < n c < Fi2- 



Appendix C. Deviation from Fraction 

To generalize the discussion in Sec. HI let us consider a 
minimum of v(a) around a mediant 



a 



m + p 



(C.l) 



of reduced fractions m/n and p/q i | Appendix A[ ). As in 
Eq. (1171) . the minimal model for this purpose is given by 

v(a) — v n V(na — m) + v q V{qa — p). (C.2) 

Differentiating this with respect to a, and substituting 



/ qm-np 

<t v i v \ : ) ■ ( c - 3 ) 



Eq. (jC.ip . we obtain 

r) ( np — mq 



v' (a) = nv n V' 



V n + q 



On physical grounds, it is natural to assume that V(x) 
is an even function, V(x) = V(—x). Then, because of 
V'(x) = -V'(-x), 



v' (a) = (qv q - nv n )V' 



, / mq — np 
V n + q 



(C.4) 



By definition, V(a) = V'(a) = for a > X. Hence, we 
get v'(a) = for X < X, where X = \mq — np\/(n + q). 
Without loss of generality, we may assume 



m p 
— < -. 

n q 



(C.5) 
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Table B.3: The main sequence with the limit divergence angle a = r 2 . 



n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Pn 


1 


1 


o 


q 
o 


C 




Q 

o 




01 




■_) 


q n 


2 


3 


5 


8 


13 


21 


34 


55 


89 


144 


n 


1 


2 


3 


5 


8 


13 


21 


34 


55 


89 


An 








1 


2 


4 


7 


12 


20 


33 


54 


Act 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


2 


6 


15 


40 


104 


273 


714 


1870 


4895 


X 


1 

2 


1 

3 


1 

5 


1 

8 


1 

13 


1 

21 


1 

34 


1 

55 


1 

89 


1 

144 



Table B.4: Fibonacci sequence. 
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1 


2 


3 


4 


5 6 


7 8 


9 


10 


11 


12 


F n 


1 


1 


2 


3 


5 8 


13 21 


34 


55 


89 


144 



Using mq — np = 1 (|9(), we obtain 

1 



X = 



n + q 



(C.6) 



for (m + p)/(n + 5). TTie optimal width X for a fraction 
p/q is given by X = l/q. For X < X, the potential v(a) 
is nearly constant for 



m X P X 

1 < a < , 

n n q q 



the width of which is 



! - ) An. 

X 



(C.7) 



(C.8) 



where we used Eqs. (|A.2[) and (IC6|) . The optimal width 
X = l/q n for the main sequence is given in Table IB. 31 
Strictly speaking, the minimum is not at a = a in 



Eq. (IC.ip . but lies at a — a + 5a by which a small cor- 
rection 5a is defined. Let us find 5a by the condition 
v'(a + 5a) — 0. In the lowest order in 5a, we have 

v' (a + 5a) = (qv q -nv n )V (X) + (n 2 v n + q 2 v q )5aV" (X) . 

(C.9) 

Hence we find 



5a = 



Vn ~ qVq I -V (X) 

n 2 v n + q 2 v q \ V" (X) 



(CIO) 



By the assumption that V(x) is localized with a half width 
X around x = 0, it is physically reasonable to suppose 
V (X) < and V" (X) > at a tail of V(x). Accordingly, 
the sign of 5a is determined by the difference between the 
denominators of the parent fractions m/n and p/q. Un- 
der a weak condition that v n decreases more rapidly than 
1/n, around a mediant a = [m + p)/(n + q) between par- 
ent fractions m/n and p/q, a minimum of v(a) is slightly 




Figure C.21: A 'vascular leaf-trace primordium' of an 8-parastichy 
(dashed line) for a = 3/8 + 5a (8a = 0.003, n c = 6) is shown along 
with a parastichy pair (3,5) (solid lines). The sign of 8a coincides 
with the sign of 2/5 — 3/8 >0, thereby the divergence angle a is 
slightly shifted from the PF 3/8 toward a limit a = 0.382. 



shifted from a — a to the side of m/n or p/q with a larger 
denominator. This is consistent with Levitov's maximal 
denominator principle on a Farey tree ([3|)- For instance, 
the position of the minimum in Fig. [16] is slightly shifted 
from 3/8 to the side of 2/5 (instead of 1/3). The effect of 
5a is neglected in Sec. [3] It is surely negligible mostly as 
confirmed explicitly. Still there are cases where the devi- 
ation 5a has discernible effects. First of all, non-zero 5a 
deforms an 'orthostichy' a = (m+p)/(n + q) into a paras- 
tichy ( | Appendix A[ ). According to the above theorem, the 
sign of 5a = a — (m + p)/(n + q) is given by the sign of 
m/n— (m+p)/(n + q) (n > q) (Fig. IC.2~Tj) . The deviation 
from a fraction 5a tends toward a limit divergence angle 
(Sec. Si]). 

Dynamical models aim to derive a = 1/t 2 determin- 
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istically ([1 H H, 0). Their results are consistent with 
each other. The original idea of the dynamical mecha- 



nism can be traced back to [26| and [22|. Most of them 



are based on a geometric assumption that the interaction 
^(|?"m ~ r n\) depends only on the distance \r m — r n \ be- 
tween the two points n and m. Generally, in an anisotropic 
system like plants, the assumption does not hold true. In 
short, a leaf is not a sphere on a cylinder surface (Fig. [T|). 
If we had adopted this strong assumption, we would have 
failed to reproduce PFs. This is because 5a enforced by 
the assumption is relatively large and not negligible. To 
compare with our model, however, it must be remembered 
that dynamical models mostly follow a long tradition of 
reseach focused on shoot apices. Unlike ours, they are not 
aimed at PFs of mature leaves. 

Finally, as a possible generalization, let us mention a 
dynamical variant of our model. In the dynamical model, 
any slight deviation 5a can be effectively important to 
determine a dynamical path of a. In the main text, we 
assumed that ao is a preset constant independent of n c . In 
contrast, one may consider a model in which ao is variable 
such that a = a m - m to minimize v{a) at a growth step n c 
determines the initial value ao of the next step n c + 1, i.e., 

a a (n c + 1) = a min (rt c ), (C.ll) 

or let ao evolve continuously or adaptively along a branch 
in the tree of Fig. 0] To begin with, the initial condi- 
tion must be set for ao(l). To avoid confusion, ao of this 
model should rather be written as a. Consider what hap- 
pens when a local maximum begins to appear around a 
minimum at a = a m ; n , when n c reaches hq + An for a m i n . 
The maximum occurs just at a rational number near but 
not at a m i n . As discussed above, a m - m is shifted by a finite 
amount 5a from the position of the maximum. Therefore, 
the sign of 5a uniquely determines the next minimum to 
be chosen. Between two newborn daughter fractions, the 
one with a larger denominator is always chosen (a daugh- 
ter fraction with a larger denominator lies on the same side 
of a parent fraction with a larger denominator). Thus, we 
obtain a > 1/3 at n c = 2, then a < 2/5 at n c = 3, and 
a > 3/8 at n c = 5, then a < 5/13 at n c = 8, and so 
forth. By following the main branch of Fig. [H we reach 
a = 1/t 2 finally in the limit n c — > oo. The bifurcation rule 
of this dynamical mechanism accords with (H) in Sec. 13.31 
Nonetheless, the causal relationship between PFs and the 
golden angle is reversed here. In the dynamical mecha- 
nism, the golden angle is the effect, or the end result a(oo) 
to be obtained generally. In the static mechanism, it is the 
cause a(0) to be selected specifically. 
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